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Jaynes-Cummings-Hubbard arrays provide unique opportunities for quantum emulation as they 
exhibit convenient state preparation and measurement, and in-situ tuning of parameters. We show 
how to reahse strongly correlated states of light in Jaynes-Cummings-Hubbard arrays under the 
introduction of an effective magnetic field. The effective field is realised by dynamic tuning of 
the cavity resonances. We demonstrate the existence of Fractional Quantum Hall states by com- 
puting topological invariants, phase transitions between topologically distinct states, and Laughlin 
wavefunction overlap. 
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Quantum systems with highly correlated states exhibit 
an exponential growth of Hilbert space with the num- 
ber of particles, making the study of arbitrary states of 
even modest systems computationally intractable. This 
problem has motivated efforts in the field of quantum 
emulation [T. A quantum emulator is designed to repli- 
cate the physics of some target system. Such emulators 
require scalable and convenient state preparation and 
measurement, and control over single and many-body 
interactions. Proposals for emulation platforms include 
ultra cold- atoms, superconductors, and sup er fluids [2]- 
|4]. Another interesting platform is coupled atom-cavity 
systems and there has been significant progress to- 
wards realising this goal p!QUT3] . Here we explore the 
physics of the Fractional Quantum Hall Effect (FQHE) 
as it relates to atom-cavity systems. 

Thirty years after their discovery, the Integer [M] and 
Fractional [15 Quantum Hall Effects are still the focus 
of intense theoretical and experimental attention [16} [TT] . 
The FQHE relies on the presence of particle-particle in- 
teractions to form highly correlated states. These states 
can exhibit anyonic, and sometimes non-abelian, exci- 
tations which are explicitly non-local. As such, the in- 
vestigation of large systems suffers strongly from the ex- 
ponential explosion in Hilbert space. While there ex- 
ist exact solutions for some FQHE systems, such as the 
Laughlin ansatz[18j, these have yet to be observed di- 
rectly in experiment. For this reason, emulation of the 
FQHE, particularly emulating the strong magnetic fields 
required, has become a major topic of interest in the sci- 
entific community [19,-^ . 

Here we show the existence of FQHE states in the 
Jaynes-Cummings-Hubbard (JCH) model [SHZ] in the 
presence of an artificial magnetic field. These states 
constitute new, strongly correlated states of light. A 
JCH lattice consists of an array of coupled photonic 
cavities, with each cavity mode coupled to a two-level 
atom [see Figs. 1(a) and (b)]. JCH systems promise 
unparallelled control and readout of the full quantum 
mechanical wavefunction. The JCH model is predicted 




FIG. 1. (a) Schematic of a square JCH lattice with a constant 
effective magnetic field. Photons moving around a plaquette 
acquire a phase A(/). (b) A single mode photonic cavity with 
frequency uj coupled to a two level atom with strength f3. (c) 
Scheme for breaking TRS in photonic cavities: a potential 
V = [V^^ + V^^ sin {u'h + A(/)^)] x {x and y in units of the 
lattice spacing) is applied to the cavities (indicated by green 
arrows) by dynamically tuning uj. The phase offset, A(/), along 
y results in the synthetic magnetic field seen in (a). 

to exhibit a number of solid state phenomena, including 
Mott /superfluid phases^, semi-conductor phvsics[23j. 
Josephson effect [24], metamaterials properties [25], and 
Bose-glass phases [26], 

We begin by introducing the JCH model, and dis- 
cuss how an artificial magnetic field can be created 
in a photonic cavity system. To demonstrate FQHE 
physics, we compute the groundstates for small systems. 
These groundstates are compared to a modified Laughlin 
ansatz, and their topology investigated. 

Each cavity in the JCH lattice is described by the 
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FIG. 2. Single particle spectrum of the JCH lattice with (a) 
A = and (b) A = 3. In each case k = 1 and (3—1. The 
spectra comprise of two transformed Hofstadter butterflies. 
Color indicates the projection into the photonic modes. 

Jaynes-Cummings (JC) Hamiltonian 

H^^ =ujL^ Acr+cr- + (3{a^a + a'a"^), (1) 

where a is the photonic annihilation operator, are 
the atomic raising and lowering operators, A the atom- 
photon detuning, P the coupling energy and h = 1. The 
states |^(e),n), where n is the number of photons, and 
g{e) are the ground (excited) state of the atom, form the 
single cavity basis. commutes with the total ex- 

citation number operator, L = a^a + <j+<j~. Therefore 
the total excitations in the cavity, ^, is a good quantum 
number. The eigenstates of Eq. ([T]) are termed polari- 
tons, superpositions of atomic and photonic excitations, 
and are a function of £ and A//3. 

The JCH model describes a tight-binding JC lattice: 

N 

jjJCH = H^c^j,^J2Hr-Y: ^^Aa, (2) 

i <'i,j> 

where Kij is the tunneling rate between cavities i and j 
and the sum over (z, j) is between nearest neighbors only. 

For large detuning (|A| eigenstates separate out 

into either atomic or photonic modes. In this limit, the 
photonic or atomic mode can be adiabatically eliminated. 
Eliminating the atomic modes, the photonic mode has a 
weak Kerr- type photon-photon repulsion [8 and the ex- 
change of energy between atomic and photonic modes 
is strongly suppressed. However, virtual processes lead 
to effective interactions in the photonic and atomic sub- 
manifolds. Photons have an atomic mediated non-linear 
onsite repulsion, making the JCH model equivalent to the 
Bose-Hubbard (BH) model [27\. Atomic modes are cou- 
pled with the effective hopping rate hz^j-^ = i<iijj3'^ / /S?^^. 
As the atomic modes are restricted to two levels, this is 
effectively a hardcore boson field for atomic states, in 
contrast to the weakly-interacting photon field. 

The QHE occurs in a 2DEG in the presence of a 
perpendicular magnetic field, which breaks time rever- 
sal symmetry (TRS). Any mechanism which breaks TRS 
manifests in the Hamiltonian as a vector potential. An 
artificial magnetic field may then be realized via the in- 
troduction of some TRS breaking interaction. Constant 



rotation or linear acceleration are the classic examples, 
leading to constant synthetic magnetic and electric fields 
respectively. Cho et a/. [29 propose a scheme for TRS 
breaking in multi-mode cavities with far detuned atoms. 
Fang et a/.[30j use magneto-optical resonators in pho- 
tonic crystals to break TRS, and Koch et a/. [311 [32] have 
recently shown TRS breaking in the context of circuit 
QED, via the introduction of a special passive coupling 
element between microwave resonator junctions. Here we 
adapt the methodology proposed by Kolovsky[21 to cav- 
ity QED where photon assisted tunneling is used to break 
TRS. An electric field with both dc and ac components is 
applied across one of the lattice axes [see Fig 1(c)]. Intro- 
ducing a phase offset, A^, in the ac field along the other 
axis results in the desired complex coupling, 27ra = A0. 
The presence of the two fields also leads to modified 
strengths of A, and /3, which can be tuned appropri- 
ately to recover the Hamiltonian in Eq. (1). The depen- 
dence of these parameters on the field is non-trivial, but 
follows the prescription in [21]. Photons do not respond 
to electric fields, however, a gradient in the cavity fre- 
quency has an equivalent effect. Recently, cavities with 
tunable resonances have been fabricated [33) l34] . This is 
achieved by the inclusion of an intra-cavity Josephson 
junction, which changes the cavity boundary conditions, 
and can be tuned via a magnetic field. Transmission line 
resonator experiments [33] have shown cj^^ can be driven 
at O(IO^) times the cavity dissipation frequency. This 
ratio provides sufficient time to observe FQH physics. 

A magnetic field is defined by a vector potential A(x), 
and is introduced into the Hamiltonian via the minimal 
substitution, p ^ p — qA{x). On a tight-binding lattice, 
a vector potential A gives rise to a complex hopping rate 
Kij KijC^'^^^^^ , where 27r6>ij = j^^ A{r) • dr. As in 
the continuum case, gauge symmetry implies that the 
only physically important parameter is the total phase, 
27ra, picked up around a closed loop, where a = ^ /^^ \s> 
the fraction of flux quanta through the loop. A constant 
magnetic field in the z direction corresponds to a constant 
a for all plaquettes on the lattice. Factors of 27r in the 
phase around a loop are physically inconsequential, so 
we only need consider a G [0, 1]. In a lattice, the effect 
of larger magnetic field strengths saturates, as distinct 
from the continuum case, where the cyclotron frequency, 
proportional to the magnetic field, has no upper bound. 

Ignoring the JC term in Eq. (1), the spectrum for the 
kinetic term, is given by solutions to Harper's equa- 
tion, resulting in the famous Hofstadter Butterfly [35], 
a fractal structure, which has q bands at a = p/g, 
(p, g) G Z. For a single excitation on the JCH lattice, 
the Shrodinger Equation is 

N 

Y^E^, = Y^H'^^^ - ^K,,{a)^^]. (3) 

i i 

Substituting an eigenvector ipf- of K with energy ki into 
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Lx X Ly 




Dim(if) 


a 


Laughlin Transition(Ac) 
overlap 


i 


4x4 


2 


512 


0.25 


0.89 


NA 


ii 


5x5 


2 


1250 


0.16 


0.99 


-1.1 


in 


6x6 


2 


2592 


0.11 


0.99 


2.5 


iv 


4x4 


3 


5472 


0.37 


0.29 


-9.1 


V 


5x5 


3 


20850 


0.24 


0.98 


-3.8 


vi 


6x6 


3 


62232 


0.17 


0.99 


NA 



TABLE I. Results for systems of size Lx x Ly sites with 
Np particles. All systems have Ci — 1/2 below the transition 
strength A = Ac. Also shown is the Hilbert space dimensions 
Dim(i7), and the Laughlin overlap. 

Eq. (|3| yields a single site Hamiltonian: 

which is transformed to Hq^ by A ^ A + A:^ . Thus the 
energies and eigenstates are: 



-(A + )/2 ± + (A - )74 



When the detuning is small (A ^ 0), there are two 
squashed Hofstadter Butterflies with a gap [Fig. 2(a)]. 
As the relative JC interaction strength j3 decreases, the 
two parts converge to recover the original butterfly. 

As previously discussed, a large detuning separates out 
the atomic and photonic states, as shown in Fig. 2. The 
spectrum is a butterfly with width ±4/^ centered around 
— /3^/A, corresponding to the photonic part, and one with 
width ±4.n^//\ centered around A + /3^/A [Fig. 2(b)]. 

The FQHE occurs for systems at sufficiently low tem- 
peratures where the flux filling factor, v = Np/N^ is 
some non-integral rational u = p/q, with Np excitations 
and N(f) total flux quanta. Here, particles lie predomi- 
nately in the lowest Landau level(LLL). When there is an 
inter-particle interaction the groundstate has long range 
off-diagonal order and an energy gap. 

We study the FQHE on a JCH lattice with periodic 
boundary conditions, to avoid edge effects, and focus on 
the u = 1/2 state, which is the most stable and accessible 
fraction for bosons, compared to the u = 1/3 found in 
the electronic case. Choosing a lattice of dimensions L^ 
and Ly^ and the number of excitations, A/p, fixes a = 
2Np/LxLy. As the state space grows quickly, we are 
constrained to the small systems in Table [ij 

The simplest FQHE states are described by the Laugh- 
lin ansatz, which is an exact solution for particles in a 
magnetic field with a contact interaction, and u = 1/q. 
Thus for a lattice, where the inter particle interaction is 
only on site, such wavefunctions can be very good ap- 
proximations to the true groundstate. The excitations of 
these states have abelian anyonic statistics. The Laugh- 
lin wavefunction with periodic boundary conditions is: 




FIG. 3. (a) Energy gaps, (b) Laughlin wavefunction overlap, 
and (c) Chern numbers in the JCH FQHE as a function of the 
detuning A. For A << 0, the effective on-site energy is very 
high. In the opposite case, where A >> 0, the effective on- 
site energy is much lower than the energy from the pressure 
and the gap is due to multiple site occupancy. Dashed lines 
indicate the transition from a fractional 1/2 state to a non- 
interacting particle state as determined by the Chern number. 
Colors correspond to configurations in Tablejl] i-black, n-dark 
blue, in-red, iv-green, i;-light blue. 

where Fcm{Z) is a function of the centre of mass Z = 
^ and freii^) depends only on the relative particle 
separations: 



Fcm{Z) 
frel{z) 



Np/q^{N^-2)/2q 
i<j ^ ^ 



where is the generalized Jacobi theta function [; 



(5) 
and 



^L{z)<xFUZ)frel{z)'[[^f 



(4) 



Oi the odd theta function. are the single particle 
states e~^^/^ of the LLL. 

For bosons (fermions) q must be even (odd) so that 
^ has the appropriate symmetry. We define a version 
of the Laughlin wavefunction for states on a JCH lattice 
by replacing the single particle Landau wavefunctions i/j^ 
with a corresponding single polariton JCH wavefunction. 
This approximation allows the computation of the over- 
lap between the explicit groundstate of H'^^^{a) and 
the Laughlin ansatz. Direct product states of single JCH 
states are not defined at points with multiple exciations, 
however, as ^{z) =0 for any Zi = Zj the issue is avoided. 

Significant overlap between our numerical results and 
the Laughlin ansatz can be a good indication of the 
FQHE. In Fig. 3(b) the overlap is shown as a function 
of A for each configuration. In the hardcore boson limit 
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(A 0), our results match those found in [37]. As the 
onsite interaction is reduced, multiple site occupancy can 
occur, and overlaps with the trial wavefunction decreases. 
To fully quantify the groundstate, two additional quan- 
tities need to be computed: (i) the topology of the state, 
and (ii) the energy gap. The Chern number quantifies 
the topology and provides an unambiguous indication of 
FQHE physics^: 




where Ox^y are the generalized periodic boundary condi- 
tions on the lattice: 

4(L^)^ = ^e'^- ^y{Ly)^ = ^e'^y (7) 

with t\. y the X and y magnetic translation operators on 
particle i . Varying 9x^y induces an electric field on the 
surface of the lattice, leading to the relationship between 
the Hall conductance and Chern number via the Kubo 
formula: (Jh oc Ci (see [38] for a full discussion). Hence, 
the topology of the groundstate with periodic boundary 
conditions is directly related to the quantization of the 
Hall conductance in the edged geometry. This measure 
provides a means of classifying the groundstate when the 
Laughlin ansatz is no longer a good representation. This 
occurs when lattice effects become significant, such as 
with large a, when the single particle states deviate sig- 
nificantly from the continuum LLLs. Hafezi et a/.[37| 
have found Chern numbers for states in the BH model. 

The degenerate groundstates, ^0(^)5 define a principal 
fiber bundle over the manifold. The Chern number 
classifies the homotopy class of the fiber bundle, which 
is a topological invariant. That is, the Chern number is 
insensitive to small perturbations relative to the energy 
gap. Only if the gap closes can the transition to topo- 
logically different states occur. Explicit computation of 
Eq.([6| is computational intensive. We instead use the 
method first proposed in [39], and used in [37], which 
allows for efficient computation of the Chern number in 
the presence of degeneracies. In this method, a phase is 
defined for the ground state at each Ox,y with respect to 
two reference states. The Chern number is given as the 
signed sum of the vortices which occur at the zeros of the 
overlap with one of the reference states. 

Figures 3(a) and (c) plot the energy gap and Chern 
number respectivly as a function of A. Both indicate 
that for some lattice configurations a transition occurs 
from the FQHE state to some uncorrelated states. In 
the JCH model (and BH), the discrete lattice gives rise to 
pressure in the system. Competition between this pres- 
sure and the on-site repulsion leads to a topological phase 
transition. In the limit of weak interactions the pressure 
dominates, and the groundstate is defined by single par- 
ticle behaviour. Between these two limits the energy gap 
closes at a single point in momentum space, marking the 



transition to a fractional state. Shown in Table [J for 
a range of lattice configurations {i — vi) is the Laugh- 
lin wavefunction overlap and the location of the value of 
A at which the transition occurs. We find that in the 
case of configurations i and vi^ no such transition occurs 
due to an exact degeneracy in the single particle energy 
spectrum. Before the transition, the ground state Chern 
number is 1/2. The Chern number changes discretely 
when the gap closes, to the non-interacting state. The 
Chern number for the system is then the sum of Chern 
numbers for single particles, given by solutions to the 
Diophantine equation Ci = sq/p — {s^Ci} G Nj4Q]. 
We find this to be the case for configurations ii^vii and 
V. For in, Ci is undefined due to the presence of level 
crossings, as we have also observed in the case of the BH 
model. The question of whether these transitions persist 
in the thermodynamic limit is still an open question[37j. 

The proposed system can potentially be implemented 
in any cavity QED framework. However, circuit QED 
systems are the most promising as they exhibit the 
largest atom-photon interactions, relative to cavity Q. 
With coherence times for qubits approaching lQ//s[4T]. 
and coupling strengths ^ 10^MHz[i5^, FQHE states on 
small lattices, on the order of 15 sites, could be produced. 

We have shown that many of the phenomena associ- 
ated with the FQHE can be realistically emulated in cav- 
ity QED systems. Cavity lattices allow direct inspection 
of quantum states, which offers an unprecedented win- 
dow into the physics of the QHE and topological phases. 
The cavity QED framework also allows for very broad 
control over the system's parameters and is readily ex- 
tensible to more complicated configurations. For exam- 
ple, by including a three level atom with evenly spaced 
levels, the photons can be given an effective 3-body con- 
tact interaction, for which the well studied Pfaffian states 
are solutions [43 . These states have non-abelian statis- 
tics, and are a basis for the description of the = 5/2 
Quantum Hall plateau. An implementation of the sys- 
tem considered here is an exciting prospect for the near 
future and will provide crucial insight into the physics of 
the Quantum Hall Effect. 

A.D.G. acknowledges the Australian Research Council 
for financial support (Project No. DP0880466). 
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